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ABSTRACT 

We calculate the Equivalent Widht of the Core and the centroid energy and relative 
flux of the 1st order Compton Shoulder of the iron Ka emission line from neutral 
matter. The calculations are performed with Monte Carlo simulations. We explore 
a large range of column densities for both transmitted and reflected spectra, and 
study the dependence on the iron abundance. The Compton Shoulder is now becoming 
observable in many objects thanks to the improved sensitivity and/or energy resolution 
of XMM-Newton and Chandra satellites, and the present work aims to provide a tool 
to derive informations on the geometry and element abundances of the line emitting 
matter from Compton Shoulder measurements. 

Key vifords: Line: formation - galaxies: active - X-rays: galaxies 



1 INTRODUCTION 

Iron Ka fluorescent lines emitted in neutral matter consists 
(in the matter rest frame) of a narrow core (corresponding 
to the line photons emerging unscattered from the emitting 
region) and several Compton Shoulders (CS; see Matt et 
al. 1991; George & Fabian 1991; Leahy & Creighton 1993; 
Sunyaev & Churazov 1996), corresponding to line photons 
emerging after one or more scatterings. While the higher 
order CS are expected to be very faint, the first order CS 
(hereinafter CSl) should now be observable in many objects 
thanks to the improved energy resolution and sensitivity of 
the instruments onboard Chandra and XMM-Newton (see 
Kaspi et al. 2002; Bianchi et al. 2002; see also Iwasawa et al. 
1997 for the only pre-Chandra observation of a Compton 
Shoulder) . 

In this paper we calculate, by means of Monte Carlo 
simulations, the relative amount and centroid energy of CSl, 
as well as the Equivalent Width (EW) of the unscattered line 
photons (Narrow Core, hereinafter NC), in both transmit- 
ted and reflected spectra. In the former case, we adopt a 
spherical geometry, and explore a large interval of the col- 
umn density. For the latter case we assume a plane-parallel 
geometry and calculate CSl properties as a function of the 
inclination angle for different values of the column density 
in the perpendicular direction. The dependence of line prop- 
erties on the iron abundance is also explored. 



2 CALCULATIONS 

The simulation code employed here is described in Matt et 
al. (1991; 1997), and we defer the reader to these papers 
for details. To avoid contamination from Compton scattered 
continuum photons, we injected only primary photons with 
energies larger than the iron Ka edge. 

For each simulated spectrum, we calculated three quan- 
tities: a) the EW of the NC with respect to the primary 
(unabsorbed) emission; b) the Centroid Energy (Ec) of CSl, 
defined as; 



J^"^ EN{E)dE 
if. NiE)dE 



(1) 



where i5roc=6.2436 keV is the mimimum energy of a line 
photon after one scattering, N{E) is the number of line pho- 
tons per unit energy at a given energy, and the integral does 
not of course include the NC; c) the ratio / between the total 
number of photons in the CSl (= J^' N{E)dE) to that in 
the NC. It should be noted that with these definitions, the 
CSl may actually include also line photons scattered twice 
or more times, provided that the energy loss per scattering 
is small enough (i.e. small scattering angles) that the emerg- 
ing photon has an energy still larger than -Erec- However, the 
chosen definitions are the most useful when comparing the 
calculations with real data. 

Element abundances and photoabsorption cross sec- 
tions are those given by Morrison & McCammon (1983) 
(the iron abundance in number is 3.3xl0~^; a value about 
40 percent higher is given by Anders & Grevesse 1989). A 
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larger (smaller) iron abundance results in an increase (de- 
crease) of the EW of the iron line (see Matt et al. 1997 for 
a detailed discussion), and in a decrease (increase) of / due 
to the larger (smaller) ratio of photoelectric absorption to 
Compton scattering. To quantify these effects, we have run 
the Montecarlo code (in the optically thick case only for re- 
flected spectra) also for AFe~0.5 and 2, where Ape is the 
iron abundance in units of the cosmic value. 

The primary spectrum has been assumed to be a power 
law with photon spectral index 2. Different indices will re- 
sult on somewhat different values of the EW of the NC (see 
George & Fabian 1991). Ec and / are instead independent of 
the spectral index, at least in the optically thin case (in the 
optically thick case there should be a small dependence on 
the slope due to the energy dependence of photoabsorption, 
and therefore of the depth at which the interaction occurs.) 

Finally, it must be remarked that we assumed cold elec- 
trons, i.e. we neglected both thermal motions for free elec- 
trons and their motion in the atoms and molecules, if they 
are bound. In real situations, where these effects are consid- 
ered, the minimum energy after one scattering, -Brcc, may 
be lower than the value given above. The adopted approx- 
imation will also lead to an artificially sharp backscatter- 
ing peak in the spectrum of Figs. 2 and 5 below, and this 
should be considered when fitting real data. (For a com- 
plete discussion of these effects, see Sunyaev & Churazov, 
1996.) We note that electron temperatures are expected to 
be large in the inner regions of accretion discs, expecially in 
Galactic Black Hole candidates, but then relativistic distor- 
tion are expected to be even larger. Electron temperatures 
are instead expected to be low in e.g. circumnuclear matter 
in AGN; in these cases, effects due to the motion of bound 
electrons should instead be considered (see e.g. Fig. 4a of 
Sunyaev & Churazov 1996). 



3 RESULTS 

3.1 CSl in transmitted spectra 

We first calculated the line spectra when the cold matter lies 
along the hne of sight (transmitted spectra) . We explored a 
wide range of column densities, from 10^^ to 5x10^* cm~^. 
We assumed a spherical distribution of the matter, with the 
source of the primary photons located at the centre of the 
sphere. The density of the matter is assumed to be constant. 

We explored the case of spherical distribution of matter 
for three different values of the iron abundance Ape with 
respect to the cosmic value (as for Morrison & McCammon 
1983): 0.5 (Fig. 0, dotted lines), 1 (solid lines), 2 (dashed 
lines). The EW of the NC increases with the column den- 
sity up to about 4x10^^ cm~^, in agreement with previous 
calculations (e.g. Inoue 1989), and then decreases. The de- 
pendence on the iron abundance is linear in the tc limit 
(tc = NH/i-5 X 10^* cm~'^), and less than linear when opac- 
ity effects become important (see Matt et al. 1997). The CSl 
centroid energy and / increase up to a Compton depth of 
about unity. The following decrease is due to the increas- 
ing importance of repeated scatterings. The dependence of 
both parameters on Ape is not dramatic and, for /, almost 
negligible up to rc ~1. 

In Fig. ^ the spectrum of the CSl is shown for three 
different column densities (Ape — i-)- 




10*= lO'" 
N„ (cm-') 



Figure 1. The EW of the NC (calculated with respect to the 
primary unabsorbed emission), centroid energy of CSl and / for 
a centrally illuminated spherical distribution of cold matter. The 
dotted, solid and dashed lines refer to Apf,=0.5, 1 and 2, respec- 
tively. 



3.2 CSl in reflected spectra 

We then calculated the CSl properties for an isotropically 
illuminated plane-parallel slab. The EW of the NC (with 
respect to the primary emission only), the centroid energy 
of CSl and / are shown in Fig. ^ (solid lines) as a function 
of fi, the cosine of the slab's inclination angle. We firstly 
assumed optically thick matter (vertical column density of 
2x10^^ cm~^) and cosmic abundances. The EW of the NC 
increases with fi (as well known, see Matt et al. 1991; George 
& Fabian 1991), as does /. Ec, instead, decreases with fj,. 

Is is worth noting that in reflected spectra /, for the 
same iron abundance, is lower than in trasmitted spectra, 
provided that the column density of the line-of-sight matter 
is larger than a few x 10'^'' cm~'^. This is because the average 
optical depth, with respect to the emerging surface, at which 
the first interaction occurs is lower in the former than in the 
latter case due to the isotropic illumination. 

We then explored the dependence of the line properties 
on the iron abundance. In Fig. ^ the results for AFe=0.5 
and 2 are shown (dotted and dashed lines, respectively). As 
expected, the EW is larger and / smaller for the larger iron 
abundance, while the opposite occurs for the smaller iron 
abundance. 

Finally, we calculated the line properties as a function 
of the vertical column density (Fig. ^). While the centroid 
energy is not significantly dependent on Nh, both / and 
the EW of the core increases, as expected, with Nh in the 
optically thin regime, to saturate for Nh ^ 5x10^^ cm~'^. In 
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Figure 2. The spectrum of CSl (in the cold electrons approx- 
imation, see text) in the transmitted case for column densities 
of 5x10^2 (solid Une), SxlO^a (dotted line) and 5x10^4 cm'^ 
(dashed line). Apf,=l. 

the Compton-thin regime, moreover, these parameters are 
much less dependent on /i than in the optically thick case. It 
is worth noting that in the optically thin regime the values 
of EW and / are larger than in transmitted spectra, for the 
same Nh , largely due to the different illumination (which is 
always normal to the surface in the spherical geometry, but 
not in the isotropically illuminated slab; in the latter case 
the effective column density becomes very large for grazing 
angles). 

In Fig. ^the spectrum of CSI is shown for two different 
values of n {AFe=l)- 



4 CONCLUSIONS 

We have calculated the properties of the iron Ka Nar- 
row Cores and of the first scattering Compton Shoulder in 
both transmitted and reflected spectra, refining and expand- 
ing previous works on the same subject (Matt et al. 1991; 
George & Fabian 1991; Leahy et al. 1993; Iwasawa et al. 
1997; see Sunyaev & Churazov 1996 for a very detailed physi- 
cal description of the process) . The intensity of the Compton 
Shoulder is at most 30-40% of the Narrow Core of the line, 
with a centroid energy typically of about 6.3 keV or more. 

Prior to the launch of Chandra, the Compton Shoul- 
der was observed only in the ASCA spectrum of NGC 1068 
(Iwasawa et al. 1997), with properties consistent with an 
origin in the Compton-thick torus of this source. Recently, 
thanks to the gratings onboard Chandra, two more cases 
have been reported. Bianchi et al. (2002) found a Compton 



Figure 3. The EW of the NC, centroid energy of CSl and / for an 
isotropically illuminated, optically thick {Ng = 2 X 10^^ cm"'^) 
plane— parallel slab, as a function of fi, the cosine of the inclination 
angle. The dotted, solid and dashed lines refer to yl^ e=0.5, 1 and 
2, respectively. 



shoulder in the HETG spectrum of the Circinus Galaxy. The 
value of /, i.e. about 20%, is consistent with reflection from 
the inner wall of the 4x10^'* cm~^ torus. 

Kaspi et al. (2002) measured a Compton shoulder in 
their 900 ks Chandra/WETG observation of the Seyfert 1 
galaxy NGC 3783. The flux in CSl is 14±4 percent of the 
narrow core, which in turn has an EW of 90±11 eV. The 
line is clearly due to reflecting matter, as no evidence for 
strong cold absorption is apparent. Looking at the values 
of / in Figs. |^ and ^, it is possible to conclude that the 
reflecting matter should have a column density of at least 
~ 10^^ cm~'^ (slightly lower values are allowed only if the 
iron is overabundant). The value of the EW is less directly 
usable, as the geometry of the reflector is unknown (the 
values in this paper refer to a 27r solid angle of the matter, 
subtended to an isotropic primary source). However, as it 
seems unlikely that the solid angle is much larger than 27r, 
low values of the iron abundance are disfavoured. 
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Figure 4. The EW of the NC, centroid energy of CSl and / for an Figure 5. The spectrum of CSl (in the cold electrons approxi- 

isotropically illuminated, plane-parallel slab with vertical column mation, see text) in the reflection case for |U=0.15 (solid line) and 

density of 5 x 10^^ (dot-dashed lines), 10^^ (long-dashed lines), /i=0.95 (dotted line). Apf.=l. 
5 X 10^3 (short-dashed lines), lO^" (dotted lines) and 5 x W"^* 
(solid lines) cm~^. Apf.=l. 
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